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The quintessential description of Kondo physics in equilibrium is obtained within a scaling picture 
that shows the buildup of Kondo screening at low temperature. For the non-equilibrium Kondo 
model with a voltage bias the key new feature are decoherence effects due to the current across the 
impurity. In the present paper we show how one can develop a consistent framework for studying the 
non-equilibrium Kondo model within a scaling picture of infinitesimal unitary transformations (flow 
equations). Decoherence effects appear naturally in third order of the /3-function and dominate 
the Hamiltonian flow for sufficiently large voltage bias. We work out the spin dynamics in non- 
equilibrium and compare it with finite temperature equilibrium results. In particular, we report on 
the behavior of the static spin susceptibility including leading logarithmic corrections and compare 
it with the celebrated equilibrium result as a function of temperature. 



I. INTRODUCTION 

A. Motivation 

Scaling concepts provide an invaluable tool in modern 
Theoretical Physics and are fundamentally important for 
understanding universal behavior and phase diagrams of 
quantum many-body systems. To date, most of these 
applications have been to equilibrium problems, where 
the basic idea of scaling to focus on the low-energy exci- 
tations of the system is most naturally applicable. Non- 
equilibrium problems, like systems prepared in an initial 
non-equilibrium state, or systems in a steady state far 
from equilibrium (like transport processes between reser- 
voirs), are already conceptually problematic for a scal- 
ing analysis since the focus cannot be on the low-energy 
excitations from the equilibrium ground state alone. A 
lot of theoretical work has recently been devoted to de- 
veloping extension of scaling ideas to non-equilibrium 
problems, namely the frequency-dependent renormaliza- 
tion group the real time reno rmalization group'^^, 
the Coulomb gas representatiorP^ and the flow equation 
approach!^. 

This paper provides a detailed account of the flow 
equation approach applied to a non-equilibrium steady 
state of a quantum impurity model. Concretely, we study 
the Kondo model with a dc-voltage bias that produces a 
stationary current. This non-equilibrium Kondo model 
is a natural candidate for our analysis since it i) can 
be realized experimentally in quantum dot experiments 
and ii) its equilibrium version is the paradigm model 
for strong-coupling impurity physics in condensed mat- 
ter theory. The current paper is a substantial extension 
of work previously presented in Ref.'^ In particular, we 
study in detail the spin dynamics and the behavior of 
the static spin susceptibility Xq i^i non-equilibrium. We 
extend previous results in the literature, which, for ex- 
ample, allows us to compare the leading logarithmic cor- 
rections in non-equilibrium with the equilibrium result 
for xo{T). 



The experimental motivation for our study is the ob- 
servation of the Kondo effect in the Coulo mb bloc kade 
regime of quantum dots, first realized in 199^2lIoIIIl, jf a 
quantum dot that is weakly coupled to two leads is tuned 
into the Coulomb blockade regime such that it carries a 
net spin, resonant tunneling through the dot leads to a 
Kondoesque increase of the conductance up to the uni- 
tarity limit upon lowering tem peratu re^^. This was first 
predicted theoretically in Refs.l^^Ell. The case of small 
voltage bias between the two leads, F ^ Tk, where Tk 
is the equilibrium Kondo temperature, can be analyzed 
using linear-response theory from the well-understood 
equilibrium ground stat d^^ '. 

The case of intermediate voltage bias, 1^ w Tk, that 
matches the linear-response to the large voltage bias 
regime, has until recently been out of reach for any 
controlled theoretical investigation. Very recently, new 
methods like the scattering state numerical renormaliza- 
tion groupP^, the time-dependent density renormaliza- 
tion groupl^ni^ and the scattering state Bethe ansatJ^ 
have been developed that can access this crossover 
regime, though much more work needs to be done be- 
fore a complete picture will emerge. 

In this paper we study the situation of large volt- 
age bias, V ^ Tk, where one expects to find weak- 
coupling physics and therefore the possibility to do a 
controlled renormalized perturbation expansion. Kamin- 
ski et al. first developed a scaling picture of the large 
voltage bias Kondo model based on the invariance of the 
current under the RG-floW^. Subsequently, Rosch et 
al. developed a more sophisticated approach based on 
frequency-dependent vertices and Keldysh diagrammat- 
ics^ 2 iTHlJ^ both approaches decoherence effects due 
to non-equilibrium spin relaxation processes generated 
by the stationary current play a key role: The result- 
ing decoherence rate Fj-oi c>c V/\v?{V/T\i)^\s essential 
for cutting off inter-lead scattering processes that are not 
immediately cut off by the voltage bias. The decoherence 
rate Fid is therefore responsible for actually making the 
situation y ^ Tk a weak-coupling problem. Such deco- 
herence effects related to the noise produced in station- 
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FIG. 1: Left: Conventional scaling picture where states are 
integrated out around the two Fermi surfaces with voltage 
bias V (here depicted for cutoff Arg < V). Right: Flow 
equation approach. Ifere all scattering processes with energy 
transfer |Ai?| < Af^q are retained in i/(Afeq). 



ary non-equilibrium states are expected to be generically 
important in non-equilibrium problems. 

The flow equation approach applied to the large volt- 
age bias regime permits an analysis of this regime that 
is complimentary to the other approaches mentioned 
above: Here the important decoherence effects emerge 
in a Hamiltonian scaling framework. We will improve 
the accuracy of some quantities using this approach, for 
example regarding the non-equilibrium dynamical spin 
susceptibility and the behavior of the static spin suscep- 
tibility as a function of the voltage bias. In addition, the 
Hamiltonian scaling picture developed here offers insights 
into non-equilibrium scaling in general. 

The basic idea of the flow equation approachP^'^ZHI! jg 
to make a Hamiltonian increasingly more band-diagonal 
by a suitable sequence of infinitesimal unitary transfor- 
mations. Band-diagonality is here measured in terms 
of the energy transfer of scattering processes, therefore 
in the flow equation framework the flowing Hamiltonian 
H(A{cq) at the "scale" Afeq only contains interaction ma- 
trix elements with energy transfer \AE\ < A^q. This 
should be compared with the conventional scaling ap- 
proach where one integrates out Hilbert space states with 
an energy larger than some cutoff Argi Siiid then suc- 
cessively lowers this cutoff. The essential idea behind 
both approaches is to organize a perturbative expansion 
in such a way as to first deal with large energy denomi- 
nators and to avoid small energy denominators; this se- 
quence allows for stable expansions even when dealing 
with nonperturbative energy scales like the Kondo tem- 
perature Tk. 

However, one important conceptual difference between 
these two approaches is that the flow equation approach 
retains all states in the Hilbert space but decouples them, 
while the conventional scaling approach actually elimi- 
nates high-energy states. If one only focuses on the low- 
energy physics in an equilibrium problem, this difference 
is unimportant and the flow equation approach is consis- 
tent with the conventional RG-flow with the identifica- 
tion Afcq oc Arg (ex plicit examples for this observation 
can be found in Ref.l^^lll). In a non-equilibrium situa- 
tion like the Kondo model with voltage bias the differ- 
ence is, however, more fundamental. Fig. [T| schemati- 
cally depicts the two schemes once the cutoff is reduced 
below the scale set by the voltage bias, A < V. In the 



flow equation framework the Hamiltonian ff(Afcq) still 
describes the stationary current flowing across the dot, 
while scaling around the separate Fermi surfaces of the 
two leads does not allow for energy-diagonal transport 
processes across the dot anymore. Therefore a naive 
application of the conventional scaling approach pur- 
ports to eliminate scattering processes with small energy- 
denominators \AE\ <C Arg across the dot, which vio- 
lates the spirit of renormalized perturbation theory. We 
will later see that it is precisely this point which allows 
us to understand qualitatively and quantitatively non- 
equilibrium decoherence in a Hamiltonian framework. 



B. Outline 

This paper is organized as follows. In Sect. II we de- 
fine the model and briefly review its conventional scaling 
analysis. In Sect. HI. A we first introduce the flow equa- 
tion method, and then in part III.B apply it to a gen- 
eral Kondo model Hamiltonian with an arbitrary num- 
ber of leads with or without voltage bias. This is the 
core part of our calculation. The resulting equations 
are then analyzed in equilibrium in part III.C, where 
we reproduce the conventional two-loop /3-function. In 
part III.D and III.E we then specialize to the case of 
interest in this paper; the Kondo model with two leads 
and voltage bias V. Using these flow equations, we work 
out and analyze the scaling picture of the Kondo model 
with voltage bias in Sect. IV. Sect. V discusses our re- 
sults for the spin dynamics in non-equilibrium, which 
are compared with the well-known equilibrium behavior 
at nonzero temperature. In particular Sects. V.C and 
V.D contain the key new results of this work regarding 
the spin dynamics and the static spin susceptibility in 
non-equilibrium. A summary of results and conclusions 
is presented in Sect. VI. Appendix A contains important 
commutators and normal-ordering results that are fre- 
quently needed, and Appendix B contains details of the 
numerical solution of the flow equations. 



II. THE KONDO MODEL WITH VOLTAGE 
BIAS 

A. Model 

We are investigating a quantum dot that can be mod- 
elled by a spin-1/2 degree of freedom S coupled to con- 
duction electrons in a left (1) and a right (r) lead with 
no external magnetic field 

H = (^p ~ f^a)'ipa'^apa~^ ^a'a ^ ' ^(a'p'){ap) ■ 

a,p,a a' ,a p' ,p 

(1) 

Here a', a = /,r label the two leads, are momen- 
tum labels, and the chemical potentials are given by 
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/i; r = ±F/2. The conduction band electron spin op- 
erators jire defined by S(a'p')(ap) = 5 Ea,/3 cl-p-^-r^^c^p^ 
where a are the Pauh matrices. The couplings Ja'a de- 
scribe the exchange interaction with the localized spin 
degree of freedom, with Jir = Jh for hermiticiy. If the 
quantum dot can be described by an Anderson impurity 
model with tunneling rates Ti^r from the left /right lead, 
the coupling constants of our effective Kondo model are 
related by Jf^ — JuJrr and Ju/ Jrr = ^i/^r^^- We define 

r Ti/Tr as the asymmetry parameter of our model. 
Wc will derive the scaling equations without making any 
of the above assumptions so that also more general com- 
plex quantum dot structures can be investigated (e.g., 
in double dot systems the relation Ji^. = JiiJ„ is not 
necessarily fulfilled). We assume a linear dispersion re- 
lation leading to a constant conduction band density of 
states p and introduce the dimensionless coupling con- 
stants gi = pJii,gr — Jrr and gt = pJir (t for transport). 
The band cutoff is denoted by D and we are only in- 
terested in the universal behavior on energy scales much 
smaller than D. 



B. Conventional Scaling Analysis 

Below we will sum up the results of the conventional 
scaling analysis for the non-equilibrium Kondo model in 
the spirit of Fig. [l] that is reminiscent of Anderson's 
"poor man's scaling" in equilibriunJSHl. The derivation 
of the one-loop results is straightforward, see Refs.'i^'SSl 
Upon lowering the cutoff Arq around the two Fermi sur- 
faces, but while Arq ^ V still holds, one finds the equi- 
librium scaling equations 



dgi 



din Arg 

dgr 
dlnARG 

dgt 
din Arg 
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plus third order terms in the coupling constants. Once 
the cutoff is reduced below the voltage bias, Arg < V, 
the strong-coupling scaling of the coupling gt stops since 
there is no sharp Fermi surface for transport processes 
on this scale. Similarly, only the inter-lead scattering 
processes still contribute to the scaling equations of gi 
and gr 



dgi 



din Arg 

dgr 
din Arg 

dgt 
din Arg 
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= ~9r 



^ 



(5) 
(6) 
(7) 



We now assume that the Kondo model can be derived 
from an underlying Anderson impurity model, and one 



easily solves the equation for gt down to the infrared 
limiiPl 



RG 



0) 



Ti + Tr HV/Tk) 



(8) 



where Tk = D cxp(— I /gi+gr) is the Kondo temperature 
of the equilibrium model (i.e., for V — 0). The IR- 
coupling ([s]) determines the current I across the dot and 
a conventional second order Keldysh calculation in the 
renormalized quantities yield^^^ J = (37r/4) ^(/^(Arg = 
0). This leads to perturbative result for the differential 
conductance G{V) valid in the limit V ^ Tk 



G{V) = a 



16 ln'{V/TK) 



(9) 



where 



Gu — 



2e2 

h (Ti + Tr)^ 

2e2 4 

X (l + r)(l + r-i) 



(10) 



is the conductance in the unitarity limit. 

Notice that there is no IR-cutoff mechanism for the 
inter-lead scattering processes in the scaling equations 
([5||6| for gi,gr, which therefore yield strong-coupling di- 
vergences even for large voltage bias^. However, Rosch 
et al. have pointed out the importance of spin relaxation 
processes due to the stationary current.l^ Second order 
perturbation theory in the renormalized quantities yields 
the decoherence rate 



Frcl OC 



V 



ln'{V/T^) 



(11) 



which cuts off the inter-lead strong-coupling flow and 
eliminates the possibility of two-channel Kondo physics 
in Kondo dots that can be derived from an underly- 
ing Anderson impurity model. The observation that the 
non-equilibrium Kondo model becomes a weak-coupling 
problem for Fid 3> Tk has then been exploited to de- 
rive a number of physically relevant quantities like the 
conduction electron T-matrix, the non-equilibrium mag- 
netization and the condu ctance b oth with and without 
an external magnetic field .'i'^llllisi 



III. FLOW EQUATION ANALYSIS 



A. Flow Equation Method 



The basic idea of the flow equation approachP£l232ll 
is to make a many-particle Hamiltonian increasingly 
more diagonal through a sequence of infinitesimal uni- 
tary transformations. Such a flow can be generated by 
the differential equation 



dH{B) 
dB 



[ri{B).H{B)] 



(12) 
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with some suitable antihermitean generator rj^B); 



Eq. ( 12 1 then generates a one-parameter family of uni- 



tarily equivalent Hamiltonians H{B). We set H{B = 0) 
as the initial Hamiltonian and want H(B = oo) to be the 
final diagonal Hamiltonian. In order to generate a sta- 
ble expansion it is of fundamental importance to prop- 
erly deal with energy-scale separation during the flow, 
similar to conventional scaling approaches. For small B 
(initial phase of the flow) we will decouple modes with 
large energy differences, while later for large flow parame- 
ter B we will deal with increasingly more energy-diagonal 
processes. This sequence of transformations is generi- 
cally generated by the "canonical" generator suggested 
by Wegneil^fil 



il{B)^^'[Ho{B),H,, 



-XB)] 



(13) 



where Ho{B) is the diagonal part of the Hamilto- 
nian H{B) and Hint{B) its interaction part. The flow 
parameter B then has the dimension (Energy)"^, and 

Afcq =' B^^/^ is the flow equation energy scale that ex- 
presses how energy-diagonal the Hamiltonian H{B) has 
become. Interaction matrix elements with energy trans- 
fer |Ai?| > Afcq are eliminated in 7J(Afoq), while pro- 
cesses with |Ai?| < Afoq are still retained. We will use 
the notation with B or Afcq interchangeably. 

The flow equation approach has been successfully ap- 
plied to numerous equilibriu m man y-body problems, like 
dissipative quantum system^^MlJ the two-dimensional 
Hubbard modeP^l, low-dimensional spin systemj^^ a,nd 
strong-coupling models like the sine-Gordon modeF ■ or 
the Kondo mode±^^S-2fi-. Due to its intrinsic energy scale 
separation the flow equation method generates the same 
IR-scaling flow as conventional renormalization group 
methods, with the identification A^q oc Af^g . Notice 
that in certain strong-coupling problem^^lllll the flow 
equation approach even allows a controlled systematic 
expansion in the strong-coupling phase where conven- 
tional scaling leads to a strong-coupling divergence. 

The main conceptual difference between flow equations 
and the conventional scaling approach is that in the flow 
equation approach the Hilbert space remains unchanged 
while modes become decoupled, whereas in the scaling 
approach high-energy modes are successively integrated 
out, compare Fig. |2] This allows one to evaluate cor- 
relation functions on all energy scales within the flow 
equation frame worlPSEI!^ however, for non-equilibrium 
problems this conceptual difference will turn out to be 
even more fundamental as we will see later (compare 
also Fig. [l|. Applications of the flow equation approach 
to non-equilib rium in itial state problems have been dis- 
cussed in Refs.'^^^HnZl^ the sequel, we will explain in 
detail the flow equation analysis of the stationary non- 
equilibrium problem provided by the Kondo model with 
voltage bias. 





H(Aro=A-6A) 





M H(A(,^) 

FIG. 2: a) Schematic picture of the conventional scaling ap- 
proach where states with energies [Arg — 5Arg, Arg] are in- 
tegrated out. b) Flow equation approach. Here all scattering 
processes with energy transfer |Ai5j > Afcq are eliminated in 
-ff(Afcq). The shaded areas in a) and b) schematically indicate 
nonvanishing interaction matrix elements. 



B. Kondo Model 

For an interacting many-body problem like the Kondo 
Hamiltonian ([T]) the flow equation approach generates 
higher and higher-order interaction terms in the com- 
mutators (fT2l and (13 1. Similar to a conventional scal- 



ing analysis, the flow therefore needs to be approximated 
by an expansion in a suitable expansion parameter that 
remains sufficiently small during the entire flow. In our 
problem we will use the run ning coupling constants J{B) 



as our expansion parameters^ and perform a consistent 
systematic expansion in orders of J{B). No further input 
except this systematic expansion is required to find both 
coherent and decoherence effects in our non-equilibrium 
problem in a systematic controlled expansion. Since we 
are performing unitary transformations (12) which do 



not explicitly depend on the ground state of the system 
(though implicitly as we will see below), it is actually ad- 
vantageous to first derive the differential flow equations 
for a general Kondo Hamiltonian with couplings to some 
arbitrary number of leads (with arbitrary chemical po- 
tentials). We introduce a general index t labelling the 
conduction band electrons and write 



H = 



(14) 



with Jt't — Jtt' for hermiticity. This will allow us to 
study both the equilibrium model without voltage bias 
and the non-equilibrium model on the same footing. In 
the equilibrium model there is only one lead, therefore 
t — p with electron momentum p, whereas in the Kondo 
model with voltage bias the label t is a multi-index 
labelling both lead left /right and electron momentum, 
t = (a, p) with a = l,r. 
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1. H{B) and Generator r]{B) 

During the flow equation procedure new interaction 
terms are generated that are not present in the initial 
Hamiltonian ( 14 1 . We will neglect the newly generated 
normal-ordered terms in 0{J^) and higher in H{B), 
which can only affect our scaling equations for the cou- 
plings J{B) in order and higher (since the commutator 
(12 1 always increases the order by at least one). We will 
later see that this order of the flow equation calculation is 
the minimum order necessary to obtain a controlled ex- 
pansion in the sense that the flowing coupling constants 
remain finite for V ^ Tk- Higher order calculations be- 
yond our present work can be understood as corrections 
to the results presented here. 

Within our order of the calculation, the flowing Hamil- 
tonian can be parametrized as 

H{B) = Y.',cU,„ + Y.Jt't(B)S-s^,, (15) 

t,a t',t 

where K^,^ u'u(-B) will turn out to be the only newly 
generated term in 0{J^). It obeys K^,^^,^{B) = 
-Ku'uu'{B) for hermiticity and initially K^,^ ,JB = 



0) = 0. In ([15| : . . . : denotes normal-ordering with re- 
spect to the ground state of the system without Kondo 
impurity. We use the normal-ordering procedure as in- 
troduced in Ref.l^Sl (for more details see the appendix). 
Notice that we do not need to e xpl icitly normal-order 
the term proportional to Jj,j in ( 15 1 since it automati- 
cally appears in its normal-ordered form. 

The generator of the inflnitesimal unitary transfor- 
mations is given by ( |13[ ) where the diagonal part of 
H{B) = Ha + Hint{B)is the conduction band kinetic 



energy Hq ^ J^t.a ^ 



t',t 



This leads to 



E rl^Lu ■■S-{s,,,xs^,J: (16) 

t' ,t,u' ,'U 



where 



„(i) 



vi'L'u = + ^u' - - e J Kt,t,u'u ■ (18) 

Here and in the sequel we will usually suppress the ex- 
plicit -B-dependence of all parameters and coupling con- 
stants in our notation. 

2. Terms m 0{J) and 0{.P) 
In order to work out the unitary flow of the Hamil- 



[rj{B),H{By\. This is the main calculational problem of 
the flow equation approach. For the current problem 
we will now proceed order by order in the coupling con- 
stant J with this calculation. 

Let us flrst evaluate [r]{B), Hq]. This is straightforward 
and leads to 



(19) 



t',t 



E (^t +^«-% - V)'7S„'„ 

t' ,t,u' ,U 



Next we need to flnd [ri{B),Hint{B)]. We first identify 
its contributions in order J^, which can only come from 

C2 = [E ^t't S ■ s,,, , E Ju'u S ■ K,u] ■ (20) 



The basic commutator [S ■ Sf, f , S ■ s^iy] is worked out 
in the appendix and given in ( A6 ) . Inserting it in ( 20 1 
yields 



C2 - z E i^t'-^t)Jt'tJu'u ■■S-{s,,,xs^J: (21) 
+ E ('t' +^t- ^^v)Jt'vJvt Hv) - 1/2) S-s^,, 

t',t,v 

E/ (^t' ~^ ~ 2e„) Jj,,y J^j : Cj/^C(Q, : 

t' ,t,v,a 

plus an uninteresting constant that does not contribute 
to the further fiow. Here we have defined the expectation 
value with respect to the non-interacting ground state 



(22) 



which is just the occupation number. No summation over 
a is implied in this definition. 



Comparing coefficients in (15 1, (19) and (21 1 we find 



dB 



{e,,~e,rj^,, (23) 



the following differential ffow equations 
dJft 

+ E(^t' + '^^v)Jt'Jvt {n{v) - 1/2) 

V 

+0(J=^) 

-(% + -e* -^uf^t'tyu 
+0(J3) 



dB 



(24) 



tonian (12 1 we now need to calculate the commutator 



From (21 1 one would conclude that a new term 

EE^*'*(^) ^^^'o^*- 



(25) 
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is generated in the Hamiltonian H{B) that has not been 
included so far 



4. Flow Equations 



dVft ^ 3 
dB 16 



We include the term (30) in the differential equation 



_l_ g _ 2£ ) J J o{J^) (26) ( 23 1 and find the following set of flow equations that com- 



However, we will see below that this contribution van- 
ishes identically due to symmetry reasons (in fact, it only 
appears in 0(J^)). 



3. Terms m 0{J^) 

Our goal is to derive the /3-function of the Kondo 
model up to two-loop order within the flow equation 
framework. This implies that we need to calculate the 
term s in order that contribute to the differential equa- 
tion ( 23 1 for J^n{B) . Notice that our differential equation 
for „/„(-B) generates this term (which is not present 
for S = 0) in order J^. Also notice that we do not 
need to calculate its 0(J'^) -contribution because this can 
only feed back into the differential equation for Jj,j(_B) in 
0{J'^). We are therefore left with identifying the terms 
with structure S-s^^ in the normal-ordered commutators 



(27) 



and 



v' .v.w' .W 



(2) 



T.Ju'uS-s,,,J ■ (28) 



We can combine these two contributions into 



X S / 



(29) 



is worked 



The commutator [S ■ s^,^ , : S ■ (s^, 
out in the appendix and given in (A12). We only need 



the contribution with structure S-s^^ and inserting yields 



(30) 



2 / y X-'-u --U' ' t' 

X {n{u') (1 - n{u)) + n{u) (1 - n{u'))) . 

We have now derived all the necessary terms to write 
down the flow equations to the desired order. 



pletely determine the flow of the Hamiltonian H{B) up 
to the resp. order in J 



dB 



(31) 



+ 2^v)Jt'vJvt Hv) ~ 1/2) 

V 

u' 

X (n(u') (1 - n{u)) + n{u) (1 - n{u'))) 

- -ih'+^u'-^t-^ufK^'tyu (32) 
+0(J3) 



In order to complete our analysis we flnally need to show 



dB 



that the terms Vt't{B) from (25 1 indeed vanish to the 



required order. This is straightforward by solving (31 1 in 
linear order 



J,,{B) = {J,,{B = 0) + 0{J^)) e-^K'-*)^ 



and inserting in ( 26 1 
3 



dB 
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Y,J,,^{B^Q)XAB^Q) 



(33) 



(34) 



x(e^, + ,,_2eJe-^(K'-.)^+(^.-*)^) 
+0(J3) 

Since we assume a linear dispersion relation and initial 
coupling constants Jt't{B = 0) that are constant (at least 
in the low-energy limit) we only need to know 

^(., + .,-2.Je-^«V-")^+(^-^')^) = (35) 



which holds up to possible band edge effects that do no 
influence the universal low-energy physics for Tk <C D. 
Therefore 



dB 



(36) 



and the Vj/t-terms do not influence the flow equations 
(31 1 and (32 1 in the orders that we are interested in. In 



fact, a more detailed analysis shows that even the 0{J^)- 
terms on the rhs of (36 1 vanish; the first nonvanishing 



contributions arise in order J 



Eqs. (31) and (32) constitute the main computational 



results of this paper and will now be analyzed in various 
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settings. Notice that they apply for a general spin-1/2 
multi-lead Kondo model where we did not yet have to 
specify whether the system is equilibrium or not. The 
derivation of these equations followed from a straightfor- 
ward application of the canonical flow equation frame- 
work as an expansion in powers of the flowing coupling 
constants. Also notice that the linear terms in (31 1 and 



( 32 1 generate the canonical exponential decay of the cou- 



pling constants (33 1 that make the Hamiltonian H{B) 
increasingly diagonal in energy space (see Fig. [2]) with 
the identification Afgq = B~^^^. 



C. Equilibrium Model 

As a consistency check and in order to gain some in- 
sight into the system of flow equations we now first an- 
alyze the equilibrium Kondo model at zero temperature 
T = 0. Here the index t labelling the conduction band 
electrons in ( 14 1 only consists of the electron momentum 



and we have one exchange coupling J. We use the fol- 
lowing approximate parametrization of the running cou- 
plings as a function of B 



pJp.p{B) = g{B)e 



(37) 



where g{B) is determined from the flow in the IR-limit 
Ep, = ep = (we choose Cp = 0). This parametriza- 
tion is asymptotically correct in the IR-limit and we can 
identify g{B) with the dimensionless running coupling 
constant of the conventional scaling approach. 



Inserting the parametrization (37 1 into the flow equa- 



tion (|3l| for e^, = = 



Eq. (41) can be solved easily 



Koo,e'eiB) = -(e'-e ) g-^^^'-^ r dB' g\B') . (42) 

Jo 

From the scaling equation in second order for g{B) one 
obtains the well-known slow logarithmic increase during 
the flow, therefore the integral in (42) is dominated by 



large values of B' and one can write 

ifoo,e'e(i?) = -(e' - e ) e-^(^'-^ B [g^B) + 0{g^)) . 

(43) 

Putting everything together in ( 38 1 and performing the 



remaining integrations we obtain 

dg ^ g" g^ 0(.g^) 
dB 2B AB B 



(44) 



In terms of the scaling parameter Afoq = B ^^'^ with 
dimension energy one can alternatively write 



dg 



din A 



fcq 



(45) 



with the correct equilibrium /3-function to two-loop 
ordei|23 



(46) 



As expected for the equilibrium model it therefore makes 
no difference whether one derives the IR-behavior with 
respect to the conventional scaling parameter Arg or the 
flow equation scaling parameter Afcq (see Fig. l2l. In the 
next section we will explore how these results change in 
the non-equilibrium situation. 



dg 
dB 



de e 



-2Be 



(-2e)(n(0-l/2) 



— OO 





/U POD 
de' / de 2(e ~ e') g e'^'^'' ' 
-OO Jo 



(38) 



where we have insert the zero temperature occupation 
numbers 



nie) = ei-e ) 



(39) 



Eq. ( 38 1 is integrated starting from B — D ^ where D is 



the conduction band width (UV-cutoff). We also need 
to the flow of the ivT-terms that are initially not present 
in the Hamiltonian but generated during the flow: from 
( 32 1 we deduce 



e'e,00 



dB 

dKoo,e'e 



dB 



= (40) 



= -{e' - ef Koo,e'e (41) 



-g2e-^(''- )^ (e'-e) 



D. Non— Equilibrium Model (general parameters) 

For the non-equilibrium setting with voltage bias V 
between the left and right lead the label t in (14) is a 



multi-index t = (a,p) labelling both left/right lead (a = 
l,r) and electron momentum p. We choose the Fermi 
energy of the left lead to be +V/2 and of the right lead to 
be —V/2, therefore we have a current flowing from left to 
right. For now we will only look at zero temperature T = 
0. 



In our flow equations (31 1 and (32) we have to differ 



entiate between the various scattering processes between 
the leads: 



dB ^"P' P {a'p'){ap} 

X] X!(V + S ~ '^^q)'^(a'p')(bq)'^(bq)(ap) 
b=l,r q 

x(nb(9)-l/2) 

E E(2^?-2e^'+S-V) V^' 



(47) 



)ibq) 



b' ,h—l,r q' ,q 
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^ (^(b'q'){bq).ia'p')iap) ^ (a' p')(ap) ,{b' q'){bq)) 

X {ub'iq') (1 - Ubiq)) + Tibiq) (1 - ni,>{q'))) 
for a'a — l,r and 



'^^{a'p')(ap),(b'q')(bq) 

dB 



(48) 



-(V - ^q)J[a'p'){ap)J{b'q')(bq) 

+0{J^) 

for 6', b,a' ,a — I, r. The effect of the voltage bias enters 
here only through the different ground state expectation 
values in the left and right lead 



ni{e)^ei-e +V/2) , n,(e ) = e(-e - F/2) 



(49) 



E. Non-Equilibrium Model (jf^ — JuJrr) 



Eqs. (47 1 and (48 1 can be simplified considerably if the 



Kondo model can be derived from an underlying Ander- 



son single impurity model. As explained in Sect. II A the 
coupling constants then fulfill the relations Jf^ = JuJrr 
and Jii/J„ = rj/Fr = r. Notice that the impurity 
physics is invariant under the exchange r <-> 1/r since 
this just amounts to an exchange of the two leads. One 
can verify easily that the impurity spin only couples to 
the following linear combination of left and right lead 
fcrmion operators 



J. dcf 
Jpa 



1 



V 1 



1 



(50) 



which obey the usual anticommuation relations 

{fpaJl'p} = ^pp'^ap- The flowing Hamiltonian 
then takes the form 

HiB) = J2'pfUpc. + Y.W^)^-'p'P (51) 



+i 



p',p 



with the conduction band spin operators s 
\ J2a.g fl'a^a0fv0 defined for the /-operators. Eqs. (47) 



p p 



and ( |48[ ) simplify to the following form where the external 
leads do not appear explicitly: 



dJp'p 
dB 



~{^p' ^p) Jp'p 



(52) 



'E(V + S - 2e,) (n^iq) - 1/2) 



x(K — K 1 

V q' q^p'p p'p^q' q) 

X {n^iq') (1 - nj{q)) -f n/g) (1 - n/g'))) 
+0{J') 



^^p'p,q'q 
dB 



-(v + 

^i^q' ^ ^q)'-^p'p'-^q'q 

+0(J3) 



(53) 



Here we have the initial condition Jpip{B = 0) = Jn + J,.^ 
and the zero temperature Fermi distribution function for 
the /-operators 



n/(p) 



if'' f ) 

\J paJ pa/ 

T-;—nr{p) 



1 + r 
1 



e > ^ 



(54) 
(55) 



< 



V 
2 



In the sequel we will restrict ourselves to the analysis of 



( 52 1 and ( 53 1 , that is we only investigate a Kondo model 



that can be realized from an underlying single Anderson 
impurity model. 



IV. 



SCALING PICTURE OF THE KONDO 
MODEL WITH VOLTAGE BIAS 

A. IR— parametrization 



Eqs. (52 1 and (53 1 contain the full information about 



the Hamiltonian flow to the resp. order in the coupling 
constant. We will analyze these equations both numeri- 
cally and analytically to understand the scaling behavior 
of the non-equilibrium Kondo model. In order to get 
some first insights, we first perform an analytical anal- 
ysis that is based on some additional approximations. 
However, we will later see that these approximations are 
justified with very good accuracy by comparison with the 
exact numerical solution. 

The flow equation differential equations are connected 
with conventional scaling equations by using an approx- 
imation of the following form 



pJp'piB)^u (B), 



(56) 



q',q 



which parametrizes the couplings in terms of the dimen- 
sionlcss running coupling constants u^{B) on an aver- 
aged energy scale e— (e^, -I- £p)/2. (In fact one can 

choose any energy in [min(ep, tp,), max(ep, e^,] with ex- 
cellent accuracy.) This ansatz solves the linear part of 
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the flow equation ( 52 1 and allows us to perform the sum- 



mations over q in closed form. One finds 



dup 
~dB 



u. 
2B 



B\l + 



l + r- 



-2B(-e^ + V/2y 



(57) 



-2j de^,de^ (e^,_g2^^g-2S(V-s)= 

X / dB' Up{B')Uq{B') 
JO 

•Kfi ^.{q') {I - n j{q)) . 

In the cubic term a whole range of values q' , q contributes 
to the integral where the product of the Fermi functions 
is nonzero. Again with very good accuracy we can replace 
the couplings Uqi , Uq under the integral by an average over 
the window [— y/2, V/2\ that is responsible for transport, 



^ J-V/2 



(58) 



The integrals over de^, and de^ can then be performed in 
closed form and one arrives at the following differential 
equation 



dUp 



(59) 



( 1 ^2B{~e^-V/2f 

2B VI + r ' 



1 + r- 



,-2Bi-e^+V/2y 



MB) 
4B2 



dB' Up{B')ut{B') 



/ r + r 
V(l + r)(l + r-i) 
1 



(l + r)(l + r-i) 



IT (2^ 



-2BV^ 



V2^V er{{V2BV)) 



We first use this to derive a closed equation for Ut{B) 
from ( 58 1 by averaging over p 



dut _ u^tiB) x/tt 



dB 2B Vssy 
Ut{B) '•^ 



erf(V2Bl/) 



(60) 



4B2 



/ dB'u^,{B') 
Jo 



(l + r)(l + r-i) 
1 

"(l + r)(l + r-i) 



2e- 



-2BV' 



V2TTBVer{{V2BV)) 



Of particular importance are also the running coupling 
constants at the left and right Fermi surfaces 

MB) = u^^=v/2iB) , Ur{B) u^^_y^,{B) , (61) 



which will turn out to be relevant for the behavior of 
the quasiparticle resonances (Kondo peaks) and which 
determine the phase diagram since they correspond to 
extremal values of Up{B — oo) as a function of B (see 
Sect. IV E I . Their flow equations follow immediately from 
(|59l): 



dui 
dB 



uf{B) 



1 



-2BV' 



2B Vl + r-1 1 
Ut{B) 



(62) 



4B2 



dB' ui{B') ut{B') 



r + r 
(l + r)(l + r-i) 
1 



dUr 



(l + r)(l + r-i) 



^^-iB)( 1 



^(2e"2i3V' + ^/2^v erf{V2BV) 



-2BV^ 



dB 2B \1 

f-B 



(63) 



MB} 
4^2 



dB'ur{B')ut{B') 



r + r 



(1 



0(1 
1 



(l + r)(l + r-i) 



-2BV^ 



V2ttBV erf{V2BV)) 



The three differential equations ( 60 1 , ( 62 1 and ( 63 1 are 



one of the key results of this work and allow us to de- 
scribe and understand the scaling behavior of the non- 
equilibrium Kondo model.f^ 



B. Scaling Analysis 



We will first analyze (60 1, (62 1 and (63 1 in the initial 



phase of the flow where Afeq ^ V. As expected, all three 
differential equations coincide and one finds for a — l,r,t: 



dUg 

dAfeq 



'a. 



feq 



U 

2A 



fcq 



OK) 

Afoq 



(64) 



Here we have used the approximation dB' u^{B') = 
B (B) , which holds plus correction terms in higher or- 
der due to the slow logarithmic flow of Ua{B). The ini- 
tial behavior of the running coupling constants is there- 



fore determined by the equilibrium /3-function (46 1 as 
expected. With our initial condition u{B — 0) — gi + 
this also agrees exactly with the scaling equations ^ 
from the analysis of Kaminski et al^ (with the identifl- 
cation Arq = A^q). 

Once the flow parameter is smaller than the voltage 
bias, Afoq <C V, the scaling equations take a different 
structure. The effective transport coupling obeys 



dut 



dA 



feq 



2 (l-hr)(l + r-i) 



(65) 
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and the couplings at the left and right Fermi surfaces, 
resp., 



dui 



dA 



fcq 



dUr 



i U 



fcq 



+Ui Ut 



dA 



fcq 



A|q 
r Afcq 

V 



A^ 

^^fcq 



TT 1 

2 (l + r)(l + r-i) 



2 (l + r)(l + r-i) 



(66) 



(67) 



One notices that the strong-coupling growth of the av- 
erage transport coupUng ut stops when Afeq ~ V due to 
the energy difference of the left and right Fermi surfaces. 
On the other hand, the couplings at the left and right 
Fermi surfaces, u; and Ur, still exhibit the typical Kon- 
doesque strong-coupling behavior in quadratic order on 
the rhs of ( 66 ) and ( 67 1 . This reflects the inter-lead scat- 



tering processes from the left lead back into the left lead, 
or likewise for the right lead: such processes still see a 
sharp Fermi surface for all Afcq <V. It is this obervation 
that has led to the prediction of 2-channel Kondo physics 
in the Kondo model with voltage bias based on a 1-loop 
calculatioEpl. 

However, different from the experience in equilibrium 
models the third order terms in (66 1 and (67 1 can become 



more important than the second order terms even for 
small coupling constants: this is due to the more strongly 
growing V/ A"^^^ terms as compared to the conventional 
1/Afcq behavior in quadratic order in the infrared limit 
Afcq — > 0. Since the third order term has a positive 
sign, it counteracts the second order term and can avoid 
the strong-coupling divergence. Figs. [3| and [4| show the 
numerical solution of ( 60 1 , ( 62 1 and ( 63 1 for symmetric 



(r — 1) and asymmetric Kondo dots (r — 2). One can 
see that indeed all couplings remain finite for sufficiently 
large voltage bias which shows that there is no 2-channel 
strong-coupling divergence. 

Before proceeding with analyzing the resulting phase 
diagram in more detail, we will first get some more ana- 
lytical insights into the solution of the scaling equations. 
The solution of (|65| is straightforward 



(68) 



Wt(Afcq) 



•\/l + Trel/Afeq 



Here 



„ def AT— 2 TA 

Frci = V27ru^ V 



(l + r)(l + r-i) ' 



(69) 



where is the running coupling constant on the scale 
Afcq = V (since we can safely neglect T^-ei/V in the weak- 
coupling regime) . From the solution of (|64[) we then find 



U* = w(Afcq = V) 



(70) 




FIG. 3: Universal curves for the flow of it; (equal to Ur) and 
Ut for symmetric Kondo dots (r = 1). Results are shown for 
various ratios V/Tk labelling the curves from top to bottom. 




FIG. 4: Universal curves for the flow of ui (full lines), Ur 
(dashed lines) and ut for asymmetrically coupled Kondo dots 
with r = Ti/Fr = 2. Results are shown for various ratios 
V/Tk labelling the curves from top to bottom. Notice that 
ui and u,. coincide for V = 0, therefore only ui is shown. 



with the Kondo temperature defined from the equilib- 
rium model 



Tk = Dy/gi + gr e 



-l/(9l+Sr) 



(71) 



In order to understand the implication of the flow of ut 



on ui and it is convenient to first rewrite (66 1 and (67 1 
in the following equivalent form: 



dui 

d hi Afcq 
dUr 

din Afcq 



^ 

l + r- 

ul 



ui 



din Ut 
din Afcq 
din Ut 



1 



din Afcq 



(72) 
(73) 
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Since 



this implies 



Trcl 



din ut 



dlnAfcq 2 Afoq + Trol 



du, 



din Afoq 

dUr 



= Ul 



Ul 



1 



din A 



fcq 



2 Afcq - 

1 r\ei 

2 Afcq - 



(74) 

(75) 
(76) 



These equations show a remarkable transmutation of the 
third order terms in the running coupling constant into 
linear terms once Afeq ^ Trei- This transmutation is pos- 
sible due to the A^^-terms in (66 1 and (67 1 (or, equiv- 
alently, due to the _B~^/^-terms in (62 1 and (63) as op- 
posed to the equilibrium B^^-termsJT^ For A^q ^ Frci 
these linear order terms are negligible and ui , exhibit 
typical Kondoesque strong-coupling scaling. In the later 
phase of the flow the couplings become small again with 

1 /2 

""((Afoq), Ur(Afoq) C!C A^^^ if thcy are not already too large 
on the scale Frci- 

Let us analyze the meaning of this scale Fj-ci • We have 
already quoted the result from Kaminski et a/.^- that for 
V ^ T^^i the current / is to leading order given by (l9| 



1 



(l-Kr)(l + r-i) 



(77) 



From (|69| we can see that F^ci c>c / with a propotionality 
factor independent of V, Tk and r. It is therefore natural 
to associate Fioi with the spin decoherence rate generated 
by the shot noise proportional to the current. The third 
order terms in the scaling equations therefore describe 
spin decoherence due to the non-equilibrium current that 
cuts off the strong-coupling behavior of the inter-lead 
scattering processes. This observation in our Hamilto- 
nian scaling framework is in agreement with the work 
by Rosch et al}^ that such non-equilibrium decoher- 
ence processes eliminate the 2-channel strong-coupling 
divergence in the non-equilibrium Kondo model with 

72 _ T 7 

We will next use the flow equation analysis to study 
quantitatively the interplay between coherent (equilib- 
rium) strong-coupling processes and decoherence gener- 
ated by the non-equilibrium current in a systematic ex- 
pansion in terms of renormalized parameters. Before pro- 
ceeding along these lines, it is worthwile to point out 
that decoherence acts differently from temperature in 
the Kondo scaling equations: Decoherence and strong- 
coupling physics are in competition in ( 75 1 and ( 76 1 , 



whereas nonzero temperature T already eliminates the 
strong-coupling term itself, e.g. 



dg 
dlnA 



-T/A 



(78) 



for the equilibrium Kondo model ( 37 ) in both the conven 



tional scaling and the flow equation approach. We would 



also like to mention the important question of how one 
can understand within the flow equation framework that 
the third order terms in the scaling equations are really 
associated with spin decoherence. We will postpone the 
answer to this question to the discussion of the dynamic 
spin correlation function later in Sect. |VC| where we will 
work out this correspondence in detail. 



C. Weak— and Strong— Coupling Regime 

From the scaling equation for ui, Ur and Ut we can now 
determine the scaling behavior of the non-equilibrium 
Kondo model: that is we determine the regions in the pa- 
rameter space where all coupling constants remain small 
(weak-coupling regime) or where at least one the coupling 
constants becomes large (strong-coupling regime). Since 
only Ul and Ur can continue to grow below A^q ~ V ac- 
cording to (751 and (76) (as opposed to uj), these regions 
are determined by the behavior of ui for r = Fj/F^ > 1, 
or by the behavior of Ur for r < 1. Because of the triv- 
ial exchange symmetry r ^ r^^, we can focus on r ^ 1 
below without loss of generality. 

From (75) we can deduce an (approximate) condition 
for u;(Atoq) remaining small during the entire flow 



-"/(Frcl) ^ 1 



< 

-1^2 



(79) 



In the notation of the original Kondo Hamiltonian ([ij 
this is equivalent to the inter-lead scattering processes 
being not too large on the scale F^ei 



pJii 



< 



(80) 



We can approximately rewrite this condition in terms of 
renormalized quantities by explicitly integrating up the 
quadratic part of ( 66 1 . One flnds 



1-Fr-i 



1 



leading to 



V_ 



In 



V27r 



V 
Tk 



l+r 



In 



rrcl 
1^ 



(81) 



Iii^{V/Tk) (l + r)(l + r-i) 



(82) 



Comparison with the numerical solution below shows 



that ( 82 1 can be used as an approxiate condition for the 
weak-coupling regime for r > 2. Before proceeding with 
the numerical solution of the full differential equation for 
Ul, one can deduce some important analytical insights 
from (82): the "critical" value of the voltage bias V/Tk 
increases for increasing values of the asymmetry r. For 
larger values of r this is due to the fact that decoherence 
is proportional to the current, which is suppressed for 
asymmetric coupling to the leads: //Tk is maximum for 
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we;ik-coupling regime 



> 10 : 



strong-coupling regime 



2.5 

= r,/r 



FIG. 5: Scaling picture of the non-equilibrium Kendo model 
as a function of asymmetry r — Ti/Fr and voltage bias. 
The dashed line separates the weak-coupling regime from the 
strong-coupling regime (deduced from the full numerical solu- 
tion, see text). The dotted line is the approximate analytical 
result (|82l. 



It should be emphasized again that one expects a 
smooth crossover between weak-coupling and strong- 
coupling regimes. While this issue cannot be definitely 
answered using the weak-coupling tools here, there is no 
physical reason to expect anything else than a smooth 
crossover. Notice that for the equilibrium Kondo model 
at finite temperature we do know that there is a smooth 
crossover between weak- and strong-coupling regimes. 
In this context the question arises what physical proper- 
ties actually distinguish the strong-coupling regime from 
the weak-coupling regime. We will address this ques- 
tion in detail later when we discuss the behavior of ob- 
servables: suitable observables behave quantitatively dif- 
ferent in these regimes and make the weak-coupling vs. 
strong-coupling distinction both useful and relevant for 
the interpretation of experiments. For example in the 
strong-coupling regime the local density of states remains 
nearly pinned to its Friedel value at the more strongly 
coupled Fermi surface, and it only decays once the volt- 
age bias is well into the weak-coupling regime. 



= 1. The third term on the Ihs of (82 1 thus reflects 



the r-dependence of the unitarity limit of the conduc- 
tance Gu in (10 1. 



The above (approximate) analytical results are con- 
firmed by the exact numerical solution of the full differ- 
ential equations ( 62 1 , ( 63 1 and ( 60 1 for the running cou- 



plings ui , Ur and Ut ■ The definition of "strong-coupling 
regime" versus "weak-coupling regime" is necessarily not 
unique since we expect a smooth crossover between these 
regimes. The definition used in our analysis is that in 
the weak-coupling regime the couplings ui and remain 
smaller than 0.75 during the entire flow, whereas in the 
strong-coupling regime at least one these coupling be- 
comes larger than 0.75. Choosing a somehow different 
value than 0.75 (e.g. using 0.5 or 1.0 instead) does actu- 
ally hardly change the crossover line. The actual value 
of 0.75 is motivated by the observation that then the 
T-matrix reaches the unitarity limit at the resp. Fermi 
surface in renormalized second order perturbation the- 
ory. This indicates the breakdown of our perturbative 
expansion in the running coupling constant and there- 
fore limits the region where our method is reliable. The 
numerical results for this crossover line between strong- 
coupling and weak-coupling regime are depicted in Fig. [5] 

Fig. [5] shows that the strong-coupling regime extends 
to remarkably large values of the voltage bias for asym- 
metric Kondo dots (as we have already qualitatively de- 
duced from ( |82p) . By comparison with the flow dia- 
grams m Figs~|3] and g] we can see that while there is 
no 2-channel Kondo physics (in the sense of a 2-channel 
strong-coupling divergence), the crossover to the strong- 
coupling regime is essentially given by 1-channel Kondo 
physics for asymmetric Kondo dots: the couplings at the 
more strongly coupled Fermi surface dominate the Kondo 
impurity physics. 



D. Expansion Schemes and Higher Orders 

We have already mentioned the remarkable transmu- 
tation of a third order term in the running coupling con- 
stant into a linear term, compare e.g. (66) and (751 
for u;(Afoq). Decoherence therefore essentially makes 
the running coupling irrelevant with a canonical scaling 
dimension -1-1/2 for Afoq < Fid. This teaches an im- 
portant lesson regarding the well-established notion that 
higher order terms in the running coupling constant can- 
not change the results of lower order calculations qual- 
itatively if the running coupling constant remains suffi- 
ciently small during the flow (i.e., in the weak-coupling 
regime) . It was exactly this observation which eliminated 
the possibility of 2-channel Kondo physics due to diverg- 
ing inter- lead scattering processes. 

The technical reason for the transmutation in our non- 
equilibrium model is the existence of the dimensionful pa- 
rameter voltage bias V , which allows terms like V/K^^^ 
alongside the conventional 1/Afeq-terms in the scaling 
equations ( 65 1 , ( 66 1 and ( 67 1 . Therefore the possibility of 



such a transmutation always arises whenever we have di- 
mensionful parameter in a scaling problem, which should 
be of importance also for analyzing other problems and 
interpreting the results of lower order calculations. The 
unavoidable general conclusion is that the expansion or- 
der by order in the running coupling constant is not a 
systematic expansion anymore, unless we can show that 
no such "transmutation" occurs and changes the results. 
Notice that for example the equilibrium Kondo model 
at finite temperature contains the dimensionful param- 
eter temperature T, which in fact also leads to trans- 
mutation in its third order terms.l^ However, already 
the second order term loses its strong-coupling behavior 
at sufficiently large temperature according to ( 78 ) , and 



therefore no qualitative changes occur due to the third 



13 



order terms. 

It is worthwile to mention again the conceptual differ- 
ence between the flow equation framework and the con- 
ventional scaling analysis around the two Fermi surfaces, 
compare Fig. [l] For Arq ^ V the scaled Hamiltonian 
has no "knowledge" of the energy scale voltage bias any- 
more and therefore cannot generate the y/A|,Q-terms 
in third order (2-loop order). On the other hand, the 
flow equation Hamiltonian still contains all sufficiently 
energy diagonal scattering processes for Af^q ^ V (com- 
pare Fig. [2|, and therefore "knows" about the energy 
voltage bias because this is the energy window where 
scattering processes contributing to the current are pos- 
sible between the left and the right lead. It is exactly 
this observation which leads to the y/A|,q-terms in the 
flow equation scaling equations. 

In this context one should also address the effect of 
higher order terms in the flow equation scaling equations. 
One can verify that the most IR-singular behavior in 
fourth order is V/Aj^^ , which is therefore smaller than 
the third order terms in ( 65]|67 l in the weak-coupling 
regime. In fifth order it seems possible to have contribu- 
tions that modify the •\/27r-proportionality factor in our 
result for the relaxation rate ( 69 1 , however, they can- 
not change the structure of (75) or (76 1. More quanti- 



tative results would be a formidable task equivalent to 
a 4-loop calculation. These observations indicate that 
the third order calculation presented here should give re- 
liable quantitative results in the entire non-equilibrium 
weak-coupling regime. 



E. Numerical Solution 

In the previous sections we have obtained analytical 
insights into the competition of decoherence and strong- 
coupling behavior in the non-equilibrium Kondo model. 
In order to do this we have used various approximations 
starting from the full set of flow equations ( 52 1 and ( 53 1 . 
More accurate results, which also allow us to check the 



accuracy of the approximations leading to ( 65]|67 l, can 
be obtained by solving the full systems of differential 
equations. These solutions will also be important for ob- 
taining quantitative results for the behavior of dynamical 
qua ntit ies lat er i n Sect.jV] The numerical effort for solv- 
ing (52) and (53) scales with N'^, where N is the number 



of band states taken into account for the numerical solu- 
tion: details of the implementation of the numerics are 
contained in Appendix [B| 

Figs. |6] and [T] depict the flow of the coupling constants 
Mp(Afeq) = pJpp(Afcq) for a symmetric and an asymmet- 
ric Kondo dot in the weak-coupling regime. One observes 
the buildup of strong-coupling behavior at the left and 
right Fermi surface, followed by decoherence effects that 
lead to a decay of the coupling constants. In Fig. [7] one 
can also notice the asymmetry in the buildup of the res- 
onances at the two Fermi surfaces due to the stronger 
effect of decoherence at the more weakly coupled lead. 
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FIG. 6: Running coupling constants Mp(Afeq) at various points 
during the flow for a symmetric Kondo dot with V/Tk ~ 8.0. 
One observes the initial buildup of strong-coupling behavior 
until Afcq/TK ~ 3, followed by decreasing coupling constants 
even at the left and right Fermi surfaces due to decoherence: 
there is no 2-channel strong-coupling divergence. 
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FIG. 7: Running coupling constants Mp(Afcq) at various 
points during the flow for an asymmetric Kondo dot with 
r = Vi/Tr = 2 and V/Tk ~ 25.0. One observes a much 
stronger buildup of strong-coupling behavior at the more 
strongly coupled lead (Fermi surface at ep = 12.5rK). 



Figs. |8] and [9] shows a comparison of the flow of the 
coupling constants Mi(Afeq), Mr(Ateq) and Ut ( Afcq ) from 
the full numerical solution with the solution of ( 60 ) , ( 62 1 



and ( 63 1 . One can see that our approximations in going 
from N X N differential equations to just 3 differential 
equations were very accurate: the resp. solutions agree 
very well. This agreement is somehow obvious in the ini- 
tial phase of the flow where the voltage bias plays no role. 
Therefore one can safely start the full numerical solution 
at Afoq = 4:V and not waste computer time with larger 
values of the initial cutoff. The curves agree very well 
until one reaches the decoherence scale (maximum of the 
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FIG. 8: Comparison of the flow of the couphng constants 
ui and ut from the full numerical solution of ( 52 1 and ( 53 1 



with the effective equations ( 60 1 and ( 62 1 . Depicted here is 
a symmetric Kondo dot (therefore ui — Ur). The full lines 
refer to the full numerical solution and the dashed lines to 
the effective equations. Results are shown for various ratios 
V/Tk labelling the pairs of curves from top to bottom. 




coupling constants), which was important for determin- 
ing the phase diagram. The small absolute differences 
at the decoherence scale between the full numerical solu- 
tion and the analytical approximation get somehow am- 
plified if one looks at the relative error in the later phase 
of the flow. This will, however, anyway turn out to be 
unimportant for the evaluation of observables (which are 
essentially determined by the absolute values of the cou- 
pling constants at the decoherence scale). In general one 
can notice from Figs. [8] a nd [9] that the accuracy of our 
analytical approximations ( |60[ ), ( 62 1 and ( 63 1 becomes in- 
creasingly better if the coupling constants remain small 
during the flow, that is for larger values of the voltage 
bias These obsersations therefore justify a poste- 

riori our previous analytical analysis based on ( 60 1 , ( 62 ) 
and desl. 



V. OBSERVABLES 
A. Observables in the Flow Equation Pramework 

The evaluation of dynamical quantities like correlation 
functions within the flow equation framework is rather 
different from conventional many-body techniques. The 
key observation is that while the time evolution becomes 
trivial in the diagonal basis for B ~^ oo, the structure 
of observables can be very complicated in this basis In 
order to use the time evolution with respect to the diago- 
nal Hamiltonian H{B = oo), one first needs to transform 
an observable O with the same unitary transformations 



FIG. 9: Same as Fig. |8] for the flow of ui, Ur and ut for 
an asymmetric Kondo dot with r = Vi/Tr = 2. The curves 
for V/Tk ~ 10 should only be considered until Tx/Afoq « 
0.3 since the scaling flow leaves the weak-coupling regime for 
smaller flow parameters Afeq. 



as the Hamiltonian itself (12) 
dO{B) 



dB 



= [r,{B),0{B)] 



(83) 



where 0{B = 0) = O is the operator in the original basis. 
Once has found 0{B = oo), it becomes easy to work out 
its equilibrium correlation functions: We are for example 
interested in 



Ceq(i) = (GS|O(t)O(0)|GS) 

= (GS|e'^*Oe-*-^*0|GS) , 



(84) 



where |GS) is the interacting ground state of the system. 
We denote the unitary transformation that connects the 
initial with the final basis with U{B = oo), therefore 

H = H{B = oo) = U{B^oo)H{B^Q)U\B^oo) 
6 = 0{B^oo) = U {B = oo) 0{B = Q)U\B ^ oo) 

where we have introduced the notation that designates 
an operator in the transformed basis. One can (at least 
formally) construct U{B = oo) from (12 I in the following 
way: 



U{B) =^ Tb exp ^ dB' 77(B') 



(85) 
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Here Tg denotes _B-ordering in the same way as the usual 
time-ordering: the generator ri{Bl) with the largest B'^ is 
commuted all the way to the left, next comes the sec- 
ond largest, etc. Next one inserts identities W{B — 
oo) U{B — oo) — I everywhere in (84) and arrives at 



Ccq(t) = (GS| U\B = oo) e'"' 6 er'"' 6 U{B = oo) |GS) 

(86) 

A key observation is that U{B — oo) |GS) is just the 
ground state of the diagonal Hamiltonian H, which is 
therefore trivial to write down (depending on the system 
for example just the vacuum, the noninteracting Fermi 
sea, etc.) 



|GS^) = [/(B = oo)|GS) 



(87) 



The evaluation of equilibrium correlation functions there- 
fore reduces to 



Ceq(i) = (GS^I e^^* Oe-'^* O |GS^) 



Since it is usually straighforward to solve the Heisenberg 
equations of motion for an operator with respect to the 
diagonal Hamiltonian H, this equation has been used as 
the starting point for the evaluation of correlation func- 
tions in various models like the spin-boson modeP^lMJ^ 
the Kondo model'^^, etc. Before proceeding along the 
same lines we first need to reconsider the changes in this 
derivation for a non-equilibrium model like the Kondo 
model with voltage bias. 

The fundamental difference from the equilibrium sit- 
uation is that now we are not interested in evaluating 
observables in the ground state. Rather, the relevant 
question is the behavior of observables in the steady state 
that evolves after switching on the coupling between the 
left and right lead and then waiting long enough. This 
is also the key idea of the Keldysh approach and in this 
manner one circumvents the generally unanswered funda- 
mental question on how to construct a non-equilibrium 
steady state. We follow the same idea here by replacing 
( 84 1 with its non-equilibrium generalization 



Cncq(i) =^ lim C{t,t^) , 



(89) 



where 



dof 



C{t,t,,) = {^>,\Oit + t^)Oit^)\^,) (90) 

Here we assume that the thermodynamic limit is always 
talcen before sending the waiting time tw for the measure- 
ment to infinity in (89 1. Notice that in general Cncq(t) 
could depend on the initial state and in the sequel 
we will use an initial state that is directly related to the 
experimental situation: is the ground state of the 

system when the two leads are not coupled via the Kondo 
impurity, that is the noninteracting Fermi sea in both the 
left and right lead separately. 



We now use the same unitary transformations as going 
from (IMf to ((Sel) 



Cit,t^)= (91) 
= (^-,1 L/^(B = oo) e^-ff(t+t») Oe~*-f^(*+*-) U{B = oo) 
X C/t(B = oo) e'^*™ O e-'"'- U{B = oo) 1^-,) . 

The additional unitary transformation of the time- 
evolved operator in the diagonal basis just maps back 
to the original basis, where the initial state is specified. 
This forward-backward scheme has also been employed 
in other model^ where one is interested in real time 
evolution from some given initial state. The backward 
transformation part of this scheme can be implemented 
by solving (83) from i? = oo to B = 0: the initial con- 
dition at i? = oo is just the time-evolved operator in 
the diagonal basis. In order to evaluate (89) from this 



expression we next need to say something about the ex- 
plicit transformation behavior of the spin operator in the 
next subsection. 



B. Transformation of the Spin Operator 

The transformation of the spin operator S under the 
flow equation transformation has already been worked 
out in Ref.l^ and we only give the key steps here. One 
makes the following ansatz for the flowing observable 
(l83l: 



SiB) = h{B) S + iY\ lp'p{B) S X sp,p + 0{u') (92) 



t',t 



and derives the flow equations for the coefficients from 



(83) 



dh 
dB 

dB 



^(Cp, - ep) Jp'p-fpp' nf{p') (1 
p',p 

hi^p, 



nf{p)) 

£p) Jp'p ~ T ^ ^ ((^p' ^ ^u)Jp'ulup 



+(ep - <^u)Juplp'u){l - 2n/(u)) 



(93) 



Initially h{B = 0) = 1, 7p'p(V^ = 0) = and under the 
flow one always finds h{B = oo) — (see Refs.'^* ■^°). 
One easily identifies this scale where h{B) rapidly drops 
to by noticing that "fpip{B) is initially only generated in 
order u, therefore h{B) initially remains nearly constant 
since its differential equation only generates terms in or- 
der . This picture breaks down on the scale Afoq ~ T^^i 
where h{B) starts to decay algebraically oc B~^^^: 



i) Nonequilibrium (voltage bias V): 



V 



r,,,(^) = V27^^.^(Afeq - ^)(r+,)(i + ,-1) 
ii) Equilibrium (nonzero temperature): 

r,ei(T) = V2^u2(Afeq = T)T 



(94) 



(95) 
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Notice that these energy scales are identical with the flow 
scale ( 69 1 where the transmutation of the third order 



terms in the coupling constant into linear terms occurs. 
At this point it is easy to understand this observation: 
When h{B) becomes small, the spin operator (92 1 com- 
pletely decays into its entangled form 5* x Sp/p. In the 
flow of the Hamiltonian this amounts to transforming 
the marginally relevant interaction term Jtn S ■ Sft into 
the irrelevant K-term, Kt't.p'pS ■ {stu x Vlp)- This im- 



plies that the flow equations for ut,ui,Ur { 74 76 ) become 
weak-coupling for A^q < Fj-ei as discussed before. In the 
next section we will identify the energy scale where the 
decoherence effects become important in the Hamilto- 
nian flow with the width of the zero frequency peak of 
the correlation function, that is the physical relaxation 
rate. 



C. Dynamical Spin Response and Correlation 
Functions 

In the following we will use the shorthand 7p'p for 
lp'p{B — Qo). From (88 1 one derives the following ex- 



pression for the symmetrized spin-spin correlation func- 
tion in equilibrium at nonzero temperature: 

Coq(<) \{{SM,sm) 

p',p q',q 

x(FS|{(5'x Vp)^(S'x v,)^}|FS) 
p',p 

+nf{p')il-nf{p))) (96) 

Here |FS) is the noninteracting Fermi sea and n/(p) the 
finite temperature Fermi-Dirac distribution. 

For the symmetrized spin-spin correlation function in 
the non-equilibrium stea dy s tate ( 89 1 we need to perform 
the backtransformation (911 of Sz{tw) for large t^. Due 
to the phases exp[i(e^, — e^) t^] that the matrix elements 
jqiq acquire, one can easily see that integrating (93 1 back 



from B ~ oo to B ~ only leads to higher order con- 
tributions in the flow equations. For large times we can 
therefore say 



SM^S,{t^) (l + 0(g)) 



(97) 



This implies for the spin-spin correlation function in the 
steady state ([89| 



C'ncq(i) 



o Ip'p 



{njip) (1 - rifip')) 



p',p 



+nf{p) (1 - nf{p))) 



(98) 



plus higher order corrections in the running coupling con- 
stant. Here we now have to insert the occupation num- 



bers (55 1 for the problem with voltage bias. 



Both the nonzero temperature equilibrium result ( 96 1 



and the steady state non-equilibrium result ( 98 1 have the 



same structure and can therefore be analyzed together: 
Fourier transformation yields 



C(w) ^ J dte'"*C(i) 
p 

+nf{p + uj) (1 - nf{p))) 



(99) 



where the notation p + uj stands for the state with energy 
ep + Lo. Due to the perturbative nature of the transfor- 
mation the sum rule 



(100) 



is not fulflUed exactly. For high voltage bias or high tem- 
perature this violation is typically only of order one per- 
cent, which becomes up to ten percent for voltage bias 
or temperature of order Tk- 

The same calculation for the spin susceptibility 



x(i) = -*e(t)([5.(t),5,(o)]) 

yields the imaginary part of its Fourier transform 



(101) 



'nf{p + uj){l-nf{p))) 



(102) 



The real part follows via a Kramers-Kronig relation, in 
particular the static spin susceptibiUty xo is given by 



1 



Xo 



duj 



(103) 



Notice that the fluctuation-dissipation theorerrP^ 



X"(^)=tanh(^ 



C{u) 



(104) 



is fulfilled exactly in equilibrium. In non-equilibrium the 
fluctuation-dissipation relation is fulfilled approximately 
in the zero frequency limit if one inserts an effective tem- 
perature motivated by the comparison of ( 94 1 and ( 95 ) : 



V 



toff 



(l + r)(l + r-i) 



(105) 



Typical curves for the spin-spin correlation function 
and the imaginary part of the dynamical spin suscep- 
tibility are shown in Fig. 10 for non-equilibrium (for 



varying voltage bias) and in Fig. 11 for equilibrium (for 
varying temperature). The spin-spin correlation function 
is a symmetric function of oj with a zero frequency peak 
of width Frei. The imaginary part of the dynamical spin 
susceptibility is an antisymmetric function of uj and has 
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FIG. 10: Universal curves for the spin-spin correlation func- 
tion C(cj) and the imaginary part of the dynamical spin 
susceptibility x"i^) in non-equilibrium (symmetric coupling 
r = 1). 



FIG. 12: Universal curves for the spin-spin correlation func- 
tion C{u!) and the imaginary part of the dynamical spin sus- 
ceptibility x"i^) non-equilibrium for asymmetric coupling 
of the leads r = 2. 
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FIG. 11: Universal curves for the spin-spin correlation func- 
tion C{u!) and the imaginary part of the dynamical spin sus- 
ceptibility x"{^) ill equilibrium for various temperatures. 
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FIG. 13: Universal curves for the spin-spin correlation func- 
tion C{u!) and the imaginary part of the dynamical spin sus- 
ceptibility x"i^) ii^ non-equilibrium for asymmetric coupling 
of the leads. The voltage bias is held fixed at V = 100 Tk- 



its maximum a t x "(^ ~ ^rci)- An approximate analyti- 
cal solution of m\ yieldfiS^ (see also Ref.E^l): 



C(w) - < 



1 

Trcl 
, ,2 



M^(Afcq = \UJ\) 



for \uj\ < Trcl 

for Trci < 1^1 < T'cff 

for Teff < |c^| 



u{A, =T,V) ^ forM<r,, 
x"{co) ^ <! (107) 
'"'^ for r,, < |.| . 

Notice that the zero frequency peak of the correlation 
function is directly related to the flow scale where h(B) in 
the transformation of the spin operator starts to deviate 
noticeably from 1 since essentially all the spectral weight 
is contained in the energy interval OiT^ei) according to 
(1061. This provides the desired physical interpretation 
of Trci defined in ( 69 1 : The energy scale where the Kondo 



(106) 



coupling becomes irrelevant can be identified with the 
spin relaxation rate. 
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In Fig. 12 we show the spin-spin correlation function 
and the imaginary part of the dynamical spin suscepti- 
bility for asymmetric coupling to the leads with r = 2. 
In Fig. [13] we depict the dependence on the asymmetry 
paramter r at fixed voltage bias. Increasing the asym- 
metry parameter r leads to a smaller decoherence rate 
since the effective temperature ( 105 1 becomes smaller. In 



agreement with ( 106 1 and ( 107 1 one therefore notices that 



mainly the low frequency properties are affected with the 
maxima of the curves increasing for larger asymmetry. 



D. Static Spin Susceptibility 

One of the most celebrated results of Kondo physics 
is the behavior of the static spin susceptibility xo as 
a function of temperature. This behavior encompasses 
the screening of the impurity spin via the formation of a 
Kondo bound state upon lowering temperature. While it 
is easy to deduce the reduction of the magnetic moment 
with scaling techniques, 



Xo(T) = 



1 

4T 



1 - 



1 



HT/Tk) 



(108) 



the correct result of a finite nonzero spin susceptiblity 
at zero temperature, XaiT — 0), was first obtained in 
the seminal application of the numerical renormalization 
group by Wilson^^. Later the Bethe ansatz'*^ gave the 
complete analytical solution leading to the high temper- 
ature expansion 



Xo{T) = 



1 

4T 



1 - 



HT/Tk) 

'2(in(r/r^))2 



(109) 



O {{HT/Tk))- 



Likewise, the behavior of the static spin susceptibility as 
a function of voltage bias is of fundamental importance 
for understanding non- equilibrium Kondo physics. In the 
sequel we will focus on the case of zero temperature, zero 
magnetic field and nonzero voltage bias. The case of 
nonzero magnetic field will be discussed in a separate 
publication. The leading order result for xo{TtV) for 
large temperature or large voltage bias was worked out 
by ParcoUet and Hooley^ 



Xo{T^O,V) 



(l + r)(l + r-i) 
W 



1 



AT, 



cff 



(110) 



This naturally raises the question about the logarithmic 
corrections to this behavior, similar to the key result 
( 108 1 of equilibrium Kondo physics. This question has 



to date only been answered for the leading lii''{V/T}<i)- 
corrections, which, however, vanish if the Kondo model 
can be derived from an underlying Anderson impurity 
model (see Sect. II. A). The logarithmic terms propor- 
tional to ln(y/TK) have not yet been calculated com- 
pletely, which means that the current status of the non- 
equilibrium calculation would miss the finite temperature 
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FIG. 14: Flow equation results for the static spin susceptibil- 
ity Xo for nonzero temperature (equilibrium) and for nonzero 
voltage bias (non-equilibrium, symmetric coupling r — 1). 
The data is plotted as a function of the effective temperature 

def 

Toff = V/A in the case of voltage bias (while Teff = T for 
nonzero temperature). For comparison we show the leading 
order Bethe ansatz result (110 1 for the equilibrium suscepti- 



bility, which agrees very well with the finite temperature flow 
equation data. 



logarithmic correction in equilibrium (108). More de 



tails about these calculations and results can be found in 
j^gfgEHHI based on perturbative RG and Keldysh tech- 
niques. 

Within the fiow equation framework it is straightfor- 
ward to derive the leading term (110 1 analytically, but 



the corrections to it can only be worked out from the 
full numerical solution of the system of differential equa- 
tions derived in the previous chapters. Numerical results 
for the static spin susceptibility obtained in this way via 



( 103 1 are shown in Figs. |14|and 15 



The comparison (Fig. [T4| ) of the Bethe ansatz result 
(110 1 with the flow equation data points in equilibrium 



as a function of temperature shows very good agreement 
for large V/Tk- This demonstrates that the flow equa- 
tion calculation presented here contains the leading log- 
arithmic corrections and can address this issue also in 
non-equilibrium . 

Flow equation results for this non-equilibrium spin 



susceptibility are plotted in Fig. 14 for symmetric cou- 
pling to the leads {r = 1) as a function of the effec- 
tive temperature T^g = V/A, Eq. (1051. In this way 



a direct comparison between equilibrium nonzero tem- 
perature and non-equilibrium nonzero voltage bias is 
possible: Naively, one might expect that not only the 
zeroth order result (110) can be described by the ef- 
fective temperature (105), but that this also holds for 
the logarithmic correction in (108). This is clearly not 



the case and the logarihmic corrections to Xo(^) contain 
some genuine non-equilibrium physics.'^ Due to numer- 
ical limitations it has unfortunately not been possible to 
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FIG. 15: Flow equation results for the static spin susceptibil- 
ity xo for nonzero voltage bias for various asymmetry param- 
eters r. The asymmetry parameter increases from bottom to 
top: r = 1.0,1.4,1.8,2.2,2.6,3.0. The da ta is plotted as a 
function of T^s = V/{1 + r){l -I- r"^), Eq. |l05| . The dashed 
line is an exact result for the behavior in the Teg limit 
independent of r, see text. Lines are guide to the eye only 
and the numerical error of the datapoints is at most 10%. 
The inset shows a blowup of the small voltage bias region. 



go to larger values of the voltage bias to study the full 
crossover to the analytically known asymptotic behavior 
limv^o, Vxo{V)/A = 1/4. 

Fig. [15] contains the non-equilibrium static spin sus- 
ceptibility plotted in the same way for various values of 
the asymmetry parameter r. One can see that the log- 
arithmic corrections get even larger for increasing asym- 
metry. The values of TeffXo(?eff) clearly start to over- 
shoot the asymptotic value 1/4 for larger asymmetries. 
As opposed to the equilibrium case, the effective mag- 
netic moment in non-equilibrium first increases upon 
lowering the voltage bias. This is related to the strong 
energy-dependence of the running coupling constant in 
non-equilibrium, see also Figs. |6] and [7] The behavior 
here is consistent with what one would conjecture based 
on the analytically known In^ -corrections for a Kondo 
model that cannot be derived from an underlying Ander- 
son impurity model. In that case the energy-dependence 
of the couplings becomes even larger than here, eventu- 
ally reaching two-channel Kondo physics behavior. Like- 
wise, the leading logarith mic cor rections to the static spin 
susceptibility are positive ,'^^^1111 which is reminiscent of 
our results depicted in Fig. [15] 

The behavior of the curves for smaller values of the 
voltage bias in Fig. [l5] has to be interpreted with care, 
since the running coupling constants become larger and 
therefore the flow equation calculation less reliable. Still, 
the fact that the curves for different asymmetries con- 
verge for smaller voltage bias has a simple physical in- 
terpretation: We know that Xo{V = 0,T = 0) = w/4Tk 
with the Wilson number w — 0.413 known exactly from 
the Bethe ansatz. This result is trivially independent 
of r. This translates into a constant slope w/A in Fig. 15 
(dashed line) , which all the other curves have to approach 
in the limit V <C Ik independent of r. Notice that our 



results for the non-equilibrium static spin susceptibility 
show clear indications of this behavior. 



VI. CONCLUSIONS 

This paper has served a twofold purpose. First, we 
have worked out how the flow equation method can be 
applied to steady state non-equilibrium impurity models. 
This was illustrated in detail for the Kondo model with 
voltage bias. Thereby we have shown how Hamiltonian 
scaling ideas can be used in non-equilibrium, and how 
they help to understand the physics qualitatively and 
also quantitatively. A point of particular importance is 
the transmutation of marginal terms in the Hamiltonian 
into irrelevant terms on the energy scale set by the spin 
relaxation rate, which is generated by the steady state 
current. This is precisely the effect of decoherence within 
our Hamiltonian scaling framework. 

In addition to these methodological considerations, we 
have also derived a number of new results for the non- 
equilibrium spin dynamics in the Kondo model with volt- 
age bias. In this paper we have focussed on the spin dy- 
namics since it contains frequency resolved information 
about the non-equilibrium steady state (as opposed to 
the steady state current ([9])). For the equilibrium model 
the spin dynamics is the key to understanding Kondo 
physics, and therefore the spin dynamics is a very suit- 
able tool to illustrate the differences between equilibrium 
and non-equilibrium. In particular, we have calculated 
the static spin susceptibility in non-equilibrium includ- 
ing its logarithmic corrections for a wide range of voltage 
bias. Although we had to resort to a numerical solution 
of the flow equations for this purpose, this is an impor- 
tant step since to date an analytical calculation of this 
leading correction has not been possible (at least for the 
experimentally most relevant case of a Kondo model that 
is derived from an underlying Anderson impurity model) . 

Scaling ideas are one of the key concepts for under- 
standing equilibrium many-body systems. Likewise, we 
believe that Hamiltonian scaling ideas like the one pre- 
sented here can play an equally important role for steady 
state non-equilibrium problems. 
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APPENDIX A: COMMUTATORS AND 
NORMAL-ORDERING 



Inserting everything into ( A2 1 yields after some straight- 
forward algebra 



In this Appendix we work out in detail some key com- 
mutators that are used throughout our calculations. 

1. Commutator [S ■ s^,^ , S ■ s„,„] 

We use the following fundamental property of the 
spin-1/2 algebra 



E 



and insert it into the commutator: 



(Al) 



(A2) 



= i -.S- {s^,^ X s„,,J : 

+S-s,,J,^,{n+{t)-l/2) 
~S-s.^,,5,,^{n+{u)-ll2) 



(A6) 



Normal-ordering in the first term on the rhs of this equa- 
tion acts on the fermions only 



S ■ (Sft X s„„J 5 • : (s,,, x s^,J 



2. Commutator [: S ■ (s\,]^ x 5*2,2) ■ ^ S ■ s^,^] 



(A7) 



i 



3 
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Next we need to normal-order the fermion terms. We 
introduce the following expectation values with respect 
to the non-interacting ground state 



(A3) 



n (u) 



dcf 



which are related by n+(u) = 1 — n (it) becau se o f 
{^1, 1 } ~ 6uu' ■ No summation over a is implied in ( A3 1 , 
and the expectation values are obviously independent of 
a due to spin symmetry (no magnetic field). Then 



dcf I 



«'«-n/3 • - C;'aC„/3 " K'Jafi (") (^4) 



and : c ,,„c 



For normal-ordering of 
higher-order fermion terms we follow the prescription of 
Wcgner^^ that not only subtracts the ground state expec- 
tation values, but also expectation values with respect to 
lower-order excited states: 



dcf 



(A5) 



+ ■■ cl„c. 



For notational clarity we use labels 1', 1, 2', 2, 4', 4 
instead of v' ,v,w' ,w,u' ,u in this section. For the 
purposes of this paper we also only need to identify 
terms with the structure S ■ s^n in the above normal- 
ordered commutator. This will simplify our calculation 
considerably. 

Using (Al) it is easy to show 



S ■ {Sy^ X S2'2) '■ T S ■ 



(A8) 



1 ^ 

a=l 



(^I'l X ^2'2) ■ ' '^4'. 



3 

a,6,c— 1 



^abc^'^i'- (^I'l X ^2'2T ■ ' ^i'i} 



The first term does not contain the impurity spin and 
can therefore not contribute to the terms that we need 
to extract from the commutator. We can focus on the 
anticommutator in the second term 



3 

E 

a.b=l 



{: s\,iS2i2 ■ , ^4 



(A9) 



3 

E 

i=l 



2 ■ T ^4 



'4} 



In order to yield terms with the structure S ■ s^.^ we need 
to extract the terms with two contractions (i.e., two n^- 
terms) in ( A9 ) . Similar to ( A5 ) one shows 



c\,CiC2iC2 



C4/C4 



(AlO) 



: (0 — contraction) : -|- : (1 — contraction) 
+(5i'4(5i4' n^(l')n^(l) : Cj/Cj : 
+S2'4S24'n+{2')n~{2) : c\,c^ : 
— (5i'4(524' n+(l')n^(2) : (^2''^i ■ 
— 52'4'5i4' n+(2')n^(l) : c\,C2 ■ 
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From this expression it is straightforward to show 

3 

: Sj^,j^S2/2 : S4/4 

i=l 

= : (0 — contraction) : + : (1 — contraction) 

3 

+Sv4Si4'n+{l')n-{l) -5^,2 
+W24' n+(2')n-(2)^sj,i 
-W24' n+(l')n-(2)^4i 
-Wi4' n+(2')n-(l)^sj,2 



Combining all the like terms in ( A9 1 yields 
3 

'^abci- (^I'l ^ ^2'2)° ■ ' *4'4} 



(All) 



a,fc=l 



: (0 — contraction) : + : (1 — contraction) 
+5r4<5i4' in+{l')n-{l) + n+(l)n-(l')) 42 
-(52'4<524' {n+{2')n-{2) + n+{2)n-{2')) s'y^ 



This gives our desired result by inserting into ( A8 1 
[■ ^ ■ (^I'l ^ ^2*2) ■ ' '5' • 5*4,4] 



(A12) 



= - Sv4Si4' (n+(l')?i"(l) + n+(l)7i-(l')) 5 • 5*2,2 

S2'4S24' (7i+(2')n-(2) + n+(2)n-(2')) 5 • 5*1,1 
+normal — ordered terms with different structure 



which reduces the complexity to quadratic in the number 
of momenta and provides an excellent approximation to 
the full system. 



The formal solution of ( 53 ) is given by 



p'p,q'q 



q 

B 

X j dB e 




(Bl) 



B{e 



Using this result the flow equations for the running cou- 
pling ( 52 1 in diagonal parametrization 



Jp'p 



(B2) 



are easily rewritten to the form 

du. 



p 
dB 



^(6^-e,)u|,e-2^(^.-.) {2nj{q)-l) 
q 

- Y.^nf{q') + nf{q) - 2nf{q')nf{q)) 



-2B(e r 



(B3) 



where 



^P\P 



dB Up>{B)up{B) 



(B4) 



APPENDIX B: NUMERICAL SOLUTION 



This Appendix contains some details of the numeri- 



cal solution of the full set of flow equations ( 52 1 and 



( 53 1 . The latter are solved using a standard Runge- 
Kutta algorithnPI. 

One numerical issue arises from the fact that the K- 



couplings in ( 48 1 depend on four momenta which leads to 



a large number of differential equations to keep track of. 
While this is no fundamental problem, it would limit the 
possible resolution on a standard workstation consider- 
ably. However, one can use the following approximation 



In the remaining flow equations, for example the transfor- 
mation of the spin operator, the running coupling Jpip(B) 
is simply replaced by the rhs of (B2). The equations 



above can be easily generalized to the case of general 
parameters (47 1 and (48 1. 



Using diagonal parametrization in the numerical solu- 
tion effectively reduces the number of differential equa- 
tions to quadratic in the number of momenta N. The 
runtime is proportional to instead of iV* in the full 
set. Additionally the stiffness of the differential equations 
is reduced by removing the exponential decay from the 
flow of the running coupling. On a standard workstation 
N can be choosen of 0(1000). 
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